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A  theory  of  a  grating  stimulated  echo  (GSE)  is  developed.  The  GSE  involves  the  sequential  excita-  |  I  I  /  ^ 

tion  of  atoms  by  two  counlerpropagating  traveling  waves,  a  standing  wave,  and  a  third  traveling  wave.  /  |  J 

It  is  shown  that  the  echo  signal  is  very  sensitive  to  small  changes  in  atomic  velocity,  much  more  sensi-  LECTF  33 

tive  than  the  normal  stimulated  echo.  Use  of  the  GSE  as  a  collisional  probe  or  accelerometer  is  dis- 

cussed.  ^  AUG  3  1993  $ 


PACS  number(s):  42.50.Md.  32.80  -t 


I.  INTRODUCTION 

Echo  phenomena  are  an  effective  means  for  probing  re¬ 
laxation  processes  in  the  field-free  region  between  the  ex¬ 
citation  pulses.  In  the  usual  two-pulse  echo  [l],  the  signal 
reflects  relaxation  of  the  atomic  dipole  coherence,  while  in 
the  stimulated  echo  (SE),  it  is  atomic  state  populations 
that  can  be  monitored.  Although  the  SE  is  an  effective 
method  for  studying  decay  rates  for  atomic  state  popula¬ 
tions,  it  is  not  a  particularly  sensitive  probe  of  relaxation 
associated  with  changes  in  atomic  velocity.  Increased  sen¬ 
sitivity  to  small  velocity  changes  can  be  achieved  by  using 
the  grating  stimulated  echo  (GSE)  described  in  this  pa¬ 
per. 

The  difference  between  the  SE  and  GSE  can  be  under¬ 
stood  in  terms  of  the  Doppler  phase  diagram  for  the  GSE 
shown  in  Fig.  I.  The  ordinate  on  this  graph  represents  the 
Doppler  phase  associated  with  various  atomic  density  ma¬ 
trix  elements.  A  radiation  pulse  propagating  in  the  k 
direction  produces  an  atomic  coherence  between  ground 
and  excited  states  at  time  T \  =0.  At  a  time  T j  before  this 
atomic  coherence  has  decayed  away,  a  second  pulse  prop¬ 
agating  in  the  —  k  direction  produces  ground-  and 
excited-state  population  gratings.  The  excited-state  grat¬ 
ings  quickly  decay  away,  but  the  ground-state  gratings 
can  live  for  a  very  long  time  (limited  only  by  some 
effective  ground-state  lifetime).  Owing  to  atomic  motion, 
the  ground-state  gratings  acquire  a  Doppler  phase 
2k  v(r  —  T i)  in  this  region.  At  time  T i,  a  standing- wave 
pulse  reverses  the  slope  of  the  phase  of  the  atomic  gratings 
and  they  begin  to  rephase.  At  time  7'4  =  2Ti-  7";,  the  net 
grating  phase  is  zero;  a  pulse  propagating  in  the  —  k 
direction  applied  at  this  time  once  again  creates  an  atomic 
state  coherence.  At  time  t—T^+Ti,  the  Doppler  phase 
associated  with  the  atomic  coherence  in  the  first  and  last 
time  intervals  cancel  one  another,  and  the  GSE  echo  pulse 
appears. 

In  the  “normal"  SE,  all  pulses  are  copropagating  and 
the  standing-wave  pulse  is  absent.  As  a  consequence,  no 
population  gratings  are  produced  and  the  phase  evolution 
between  times  7";  and  T 4  would  be  represented  by  a  hor¬ 
izontal  line  in  Fig.  I . 

The  sensitivity  of  the  SE  and  GSE  to  small  velocity 
changed  Sv  in  the  time  interval  (7V7':)  is  dramatically 
different.  In  the  SE,  the  velocity  changes  modify  the  per¬ 


fect  phase  cancellation  of  the  atomic  state  coherence  only. 
This  phase  change  is  of  order  |k • ^vTTl  <  |k-  <Sv/y|.  where 
y  is  the  rate  at  which  the  atomic  coherence  decays.  On 
the  other  hand,  owing  to  the  Doppler  phase  associated 
with  the  population  gratings,  the  GSE's  sensitivity  to  ve¬ 
locity  changes  is  of  order  |k- Sv(Ti  —  7';)|  |k- ^vF;|. 

The  production  of  population  gratings  in  echo  phenomena 
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FIG.  1.  Doppler  phase  associated  with  various  density  matrix 
elements  leading  to  a  grating  stimulated  echo  (GSE).  The  four 
input  pulses  K,  (1  “  1-4)  give  rise  to  an  echo  E,  at  l~ l,.  The 
tield  propagation  vectors  are  ki  “  —  k; ”  —  k< ”k,.  There  is  a 
second  contribution  to  the  echo  signal  (not  shown)  for  which 
kj  ”ki  ”  — k.;  that  contribution  is  obtained  by  reflection  about 
the  t  axis  and  involves  density  matrix  elements  p\\.  p n.  pu,  pu. 
The  corresponding  diagram  for  the  stimulated  echo  (SF)  has 
F  1  “0  and  ki  —  k;  “ k4  “k,;  the  Doppler  phase  for  the  St;  is  con¬ 
stant  in  the  time  interval  /  j:  The  atomic  energy -level  diagram 
is  shown  in  the  inset  Each  level  decays  to  an  external  reservoir 
at  rate  yi  and  the  overall  decay  rate  of  level  3  is  E—  f'u 
+  Ei;+  yi.  The  incident  fields  drive  only  the  1-3  transition. 
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was  first  discussed  by  Mossberg  el  al.  (2], 

In  Sec.  11.  an  expression  for  the  GSE  is  derived.  The 
long-lived  ground-state  grating  exists  only  if  the  sum  of 
ground-  plus  excited-state  population  is  not  conserved 
[3,4l;  consequently,  the  atomic  model  we  choose  is  one  in 
which  the  ground  state  consists  of  two  sublevels.  In  Sec. 
Ill,  two  applications  of  the  GSE  are  discussed.  It  is 
shown  that  the  GSE  can  be  used  as  a  sensitive  probe  of 
collisions  occurring  between  ground-state  “active"  atoms 
and  perturbers.  Furthermore,  it  is  shown  that  the  GSE 
can  serve  as  the  basis  for  an  atom  accelerometer.  Such  an 
accelerometer  is  similar  to  those  proposed  by  Borde  [5] 
and  demonstrated  experimentally  by  Richie  el  al.  [b]  and 
Kasevich  and  C  hu  [7] .  An  essential  distinction  is  that  the 
f  iS I  does  not  require  a  long-lived  excited  state,  as  in  [6], 
and  is  not  connected  with  the  establishment  of  a  coher¬ 
ence  between  ground-state  sublevels,  as  in  [7], 


II.  ECHO  FORMATION 


Let  us  consider  a  gas  of  three-level  atoms  interacting 
with  a  radiation  held  E(r.f ).  The  upper  atomic  level  3  de- 
cavs  to  level  I  at  rate  To  and  to  level  2  at  rate  Ti:  (see 
Fig.  I)  The  held 


Etr.f 


ze 


j- 1 


EjS)  (l  —  Tj) 


J/'\  j  =  1,2,4 
x  (cos(A.v ),  j  =3 


+c.c. , 


(I) 


consists  of  a  number  of  pulses,  each  having  central  fre¬ 
quency  U.  Pulse  j  (7  =  1-4)  has  amplitude  E,.  and  is 
characterized  by  a  pulse  envelope  function  £,(/  “  7’,)  cen¬ 
tered  at  1  — Tj.  having  a  temporal  width  of  order  r(.  Pulse 
intervals  are  defined  by  Tlk=T,-'l\  and  T \  is  set 
equal  to  zero.  All  pulses  are  polarized  along  the  r  axis. 
The  propagation  directions  are  chosen  as  A|  =  — A:=A 
and  k  4  =  ±  k.  w  here  k  =  ll/e. 


We  wish  to  calculate  the  echo  signal  generated  at  times 
t  >  T *.  The  master  equation  for  the  density  matrix  p  (in 
the  interaction  representation)  is 


—  +  r 

dt 

p33  =  -(f/7t)[E.p]33, 

(2a) 

d  4. 

7i  +  y‘ 

P11  =r3,p33-(/77i  )[E,p]|i  +  y/H,|(r)  , 

(2b) 

d  , 

~d,+y‘ 

P22  =  T32P33+  yi  IE s(r  ) , 

(2c) 

d  4. 

h  +  y 

Pi 3  =  (i/h  )[E,p],  , , 

(2d) 

where  d/dt  —  d/dt  +  ed/dx,  r  is  the  x  component  of  the 
atomic  velocity,  V=  —  d  E(r,r)  is  the  atom-field  interac¬ 
tion  Elamiltonian,  d  is  the  atomic  dipole  moment  operator, 
r  =  r„  +  r,:+y,  is  the  upper  level  decay  rale,  y/«T  is 
some  state-independent,  effective  loss  rate  for  the  atoms, 
y~  T/2  is  the  rate  at  which  the  coherence  p\\  decays,  and 
fTjfr)  (1  =  1,2)  is  the  level  1  equilibrium  velocity  distribu¬ 
tion  in  the  absence  of  the  fields. 

Equations  (2)  are  solved  assuming  that  [A  r;  K  '• 
Ery«l.  A-MTy«l,  where  A  =  I1~  w4|  is  the  detuning  of 
the  held  frequency  relative  to  that  of  the  atomic  transition 
3  •  1  and  u  is  the  most  probable  atomic  speed.  Further¬ 
more.  the  traveling-wave  pulses  are  treated  to  lowest  order 
in  perturbation  theory.  These  restrictions  are  imposed  for 
mathematical  simplicity,  they  are  in  no  way  critical  to  the 
effect  to  be  derived  herein. 

It  is  convenient  to  carry  out  the  calculation  in  the  atom¬ 
ic  rest  frame.  The  transformation  between  the  laboratory 
and  atomic  frames  is  given  by  .\'o=.v  —  i  f.  where  ,vn  is  the 
atomic  position  at  1  =  0.  Following  a  rather  straightfor¬ 
ward  calculation  [i.8.9],  one  ariives  at  the  density  matrix 
for  f  >  74.  Keeping  only  that  part  of  the  density  matrix 
responsible  for  echo  formation  al  f=r,v  (see  Eig.  I),  one 
finds 


p(/ )  =  (//) f» ) Cl  i:/E)  exp[  —  y(  r  +  1 )  —  y/7'4s]  W|(v)04 

x  exp[/A  4(  vn+ 1/  4 )  —  i A  7  4]  Uh_0*  expl/  (A :  —  A  1 )  \  ()  —  /  (A  —  A  m ) 7  >  1  ]  -Ec.c.l 

x  cos)|di|co.s[A  ( Vo+r7'))]l  |3>(  1 1  +  H. c.  (3) 


where  0,  =  —  (</* /•.' ,!h  )J°° *,</r  g,(f )  (|d,|  is  the  area  of 
pulse  /  with  |f/,|<^l  for  /  =  1.2,4),  </=<  1 1  d  -  f  3 )  and 
r  =  /  -  7  4. 

for  an  echo  signal  havings  propagation  veettir  k, 
—  ±  Ax  =  A,x.  the  positive  frequency  part  of  the  echo  am¬ 
plitude  /-.,(/)  =2rr/A7./>,  (f )  is  proportional  to  the  positive 
frequency  part  of  the  A,  spatial  harmonic  of  the  polariza¬ 
tion  f\(i )  which,  in  turn,  is  related  to  the  atomic  density 
matrix  via 


I 

7\(f )  —nde'^’i  Ml. )  f  dxe~'k’'  f°°  <7rO|p(f  )|  I  > . 

(4) 

where  /.  is  the  length  of  the  gas  cell  along  the  \  axis  and  n 
is  the  atomic  density.  Substituting  expression  (3)  into 
(4),  setting  .v=.\o  +  r/,  and  using  the  fact  that 
cos(ucosp)  “Xm(  ~  I  )"'7:m(« )<•’""*,  where  7v(u)  is  the 
Bessel  function  of  order  q.  one  finds  the  positive  frequency 
part  of  the  polarization  is 


l\(i) 


(i7lb)m/(Ei:/l  ):expl/Ar  —  y(7':i  +  r  )  —  yi  7 \»;l(04(M7f<’ 


'r  +  0i0:0\e‘U:,r') . 


1 5a) 


where 
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P--£~  £  (-l)mi2m(!0j|)  f-Ldxo 

x  J  ^  dr  IC|  (r  )exp{/.vo(A 4  +  A 2  —  A  1  +  2mk  —  A, )  +  irl TijlA 4  +  2wA  —  A, )  +  T 43  (A 4  —  A, ) ] 

+  /r  1 7’:  1  ( A  4  +  A  2  +  2  w  A  —  A  t )  —  A5  r  ]  1  ,  (5b) 

and  P' is  obtained  from  (5b)  with  the  replacements  A 1  *  —  A  |,  A  2  -*•  —  A 2. 

Let  us  consider  expression  (5b)  for  7>.  The  spatial  phase  must  vanish  for  a  coherent  contribution  to  the  echo  field.  This 
condition,  together  with  the  fact  that  k\  =  —  A  2  =  A,  leads  to  a  value  for  the  propagation  constant.  A,  =  A4  +  2(/»  —  1  )A. 
Since  An 7"  32^1  and  r,2»7'2i  it  is  not  possible  to  get  a  nonvanishing  result  after  averaging  over  c  unless 
7-43  ~  T si! (m  —  I ).  Integral  values  w  >  2  are  allowed.  It  then  follows  that  the  phase-matching  condition  [A,  =  n/cl  can 
be  satisfied  only  for  m  =  2(£41  —  T .12).  A 4  =  —  A,  and  A,  =  A  (echo  propagation  in  the  +,v  direction).  In  this  case,  the  echo 
amplitude  is 

£.,(/)  =  -(^A7./8)wc/(r32/r)-046»20rexp[-y(7'2i  +  r)-2)'/7'32  +  /A(r-7'2i)]y4(|^il)F!(r  -7'.,  >  .  <(■) 


where  F \(i )  =  /-oodrIf'|(r  )e 

The  corresponding  result  for  the  P'  term  of  Eq.  (5a)  is 
T 43 =  £32,  A 4  =  A ,  and  A,  =  —  A.  In  this  case,  the  echo  sig¬ 
nal  propagates  in  the  —  x  direction,  having  amplitude 
£.,■(/)  “ftjlEjOJ/ftJ*  exp(2/Ar ), 


III.  DISCUSSION  OF  THE  ECHO  SIGNAL: 

APPLICATIONS 

(I )  The  structure  of  the  echo  depends  on  the  homogene¬ 
ous  width  y  and  inhomogeneous  width  Am.  Two  limiting 
cases  of  experimental  interest  are  ku»y  and  ku<y. 
The  first  condition  is  typical  for  atoms  at  room  tempera¬ 
ture  while  the  latter  typical  for  sub-Doppler,  laser  cooled 
atoms. 

For  Am»x  one  finds  that  the  GSE  has  a  temporal 
width  of  order  (Am)  -l  and  is  centered  at  £  =£4+  £21-  as¬ 
suming  £4i  =  £_i2-  The  backward-directed  echo  contains 
Ramsey  fringe  structure,  owing  to  a  factor  of  exp(2/ A£ 21 ) 
in  the  echo  amplitude,  which  is  typical  for  echoes  induced 
by  the  separated  counterpropagating  waves  (10).  Such  a 
factor  is  absent  for  the  forward  directed  echo:  however,  it 
is  possible  to  create  Ramsey  fringes  for  the  forward- 
directed  echo  by  introducing  a  dilference  57'  =  £43  —  7<2 
between  the  pulse  intervals.  If  ST  s  y~\  one  does  not 
destroy  the  echo  signal,  which  is  centered  at  f,  =£4 
+  £21  —257'.  The  forward  directed  echo  amplitude  con¬ 
tains  a  factor  of  exp(2/A(57')l  which  gives  rise  to  Ramsey 
fringe  structure. 

When  Am  y.  the  time  scale  of  the  echo  experiment 
must  be  much  larger  than  y~]  to  allow  for  inhomogene¬ 
ous  relaxation  of  the  different  atomic  velocity  subgroups. 
Consequently,  the  standard  two-pulse  excitation  scheme 
(two  pulses  separated  by  a  time  interval  of  order  y_l)  III 
does  not  lead  to  echo  formation.  Excitation  pulse  se¬ 
quences  similar  to  those  of  the  SE  or  GSE.  with 
E42»y~'.  are  needed  for  echo  formation  in  the  limit 
AM«y.  The  temporal  structure  of  the  echo  is  very 
different  than  that  for  A m  »  y.  When  A m  «  y.  the  echo 
signal  appears  as  soon  as  the  pulse  at  I  “  £4  is  applied, 
and  has  a  duration  of  order  y_l.  In  the  SE.  the  echo 
propagates  in  the  same  direction  as  both  the  last  excita¬ 
tion  pulse  and  the  free  decay  signal  following  this  pulse. 


Only  the  specific  dependence  on  the  delay  £3 2  would  per¬ 
mit  one  to  extract  the  SE  contribution  to  the  total  signal. 
In  contrast,  the  GSE  propagates  in  a  direction  opposite  to 
that  of  the  last  excitation  pulse  and  can  be  spatiallv 
separated  from  both  the  driving  pulse  and  the  free  decay 
signal.  Thus  the  GSE  offers  a  distinct  advantage  over  the 
SE  for  the  observation  of  such  a  signal. 

(2)  Aside  from  the  directionality  properties,  the  major 
advantage  of  the  GSE  over  the  SE  is  its  sensitn  ity  to 
small  velocity  changes.  The  GSE  echo  amplitude. 
£(GSE).  is  proportional  to  (see  Fig.  1)  cxpIiAH  / 

—  r )lexp[2/Ar(£43  —  £32)].  while  the  corresponding  SE 
amplitude.  £(SE).  (with  k|=lc2=k4~k  and  £:=())  is 
proportional  to  expl/'Ai- ( £;  1  -  r)l.  reflecting  the  fact  that 
the  population  part  of  the  phase  diagram  (Fig.  I)  is  hor¬ 
izontal  between  £ 2  and  £4  for  the  SE.  Any  small  change 
Sr  in  r  results  in  a  phase  shift  of  order  2A5i  ( £21  +  7 k ) 
for  the  GSE.  but  only  kSrT:\  for  the  SE  Since 
7 32 »  £21  -  the  sensitivity  to  small  Sr  is  much  larger  for 
the  GSE  than  the  SE.  This  feature  has  been  exploited  by 
Kasevieh  and  Chu  in  constructing  an  atom  interferometer 
based  or.  the  Raman  effect  (7). 

The  sensitivity  to  small  Sr  makes  the  GSE  an  ideal  tool 
for  probing  atomic  collisions.  Consider  the  limiting  case 
when  kSrT 2\  «  I  and  E 1 3 £  12 «  I  -  where  Ej  <  is  ihe  col¬ 
lision  rate  associated  with  the  atomic  coherence  pi 3.  As¬ 
suming,  moreover,  that  ££32 I,  T;  «  E  (E,‘  is  the  col¬ 
lision  rate  for  atoms  in  slate  ;  =1.2.3)  and  that  the  col¬ 
lision  kernel  IE,(r'.r)  (/  =  1.2,3),  defined  as  the  probabili¬ 
ty  density  per  unit  time  that  an  atom  in  state  1  changes 
its  velocity  from  r'  to  r,  can  be  written  as  Ujtr’.r) 

—  E;  /,  ( r  —  1  ' )  [111.  u  here  /,  (t  )  =  /,(— f )  and 

—  1 .  one  can  show  that  the  SE  amplitude  is  given  by  [  I  3) 

£  (SE )  oc  .-I  (SE )  exp(  —  2  y,  7  ) 

xexp(-Ei/A  :(Sr  £:i  *  1  .  (7) 

where 

r  ■*  1 '  a  *  7  ■*  * 

.  f  (SE) - -  +  ~ - (E,,<5r-’>,-Eiu'ir:\)  . 

i  r-  2 

(Sr:),  and  7' “£32  “£43.  The  correspond¬ 

ing  expression  for  the  GSE.  is 
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£-(GSE)oc^(SE)l(r,;/r)  +  (rJ,/r;!)[.4j(2A7)-.l,(2*r)]|exp(  — 2y,7')expi-2nT[l  -//,(D)1  , 


where  ,4,(v)  “T/tl  —  J/,(e)exp(  —  hDr/#]  and  //,(/) 
=  t  >fbdl'fdejj(()exp(2ikel').  Note  that,  if  T]  =n'. 
/i(f)“/j(f),  and  r32  “0,  population  is  conserved  for 
each  velocity  subgroup  (except  for  an  overall  decay  rate) 
and  the  SE  and  GSE  amplitudes  vanish  [4], 

In  comparing  the  SE  and  GSE  signals,  one  notes  the 
following  differences:  ( 1 )  The  SE  amplitude  depends  only 
(<Sr '),  and  not  on  the  specific  form  of  the  kernel,  while  the 
GSE  amplitude  could,  in  principle,  be  used  to  extract  the 
form  of  / |(e).  (2)  The  minimum  velocity  change  in  a  sin¬ 
gle  collision  detectable  using  the  GSE  is  of  order 
(2kD _l  >  yi/2k,  whereas  it  is  (2T/T:i)»l  times 
larger  for  the  SE.  (3)  The  minimum  collision  rate  T| 
detectable  using  the  GSE  is  of  order  (2D  “ 1  >  y;/2, 
whereas  it  2(A:  2(<5r :>i D| 2) -,3>  I  times  larger  for  the 
SE.  Thus,  the  GSE  is  a  much  more  sensitive  probe  of  col- 
lisional  effects  than  the  SE. 

For  thermal  samples,  one  can  take  y/  — 104  s_l  to 
simulate  some  effective  lifetime  of  the  atoms  in  the  atom- 
field  interaction  region.  In  that  case,  one  can  use  the  GSE 
to  measure  dr's  as  small  as  0.1  cm/s  (taking  A=I05 
cm-1)  and  collision  rates  Y\/2n  as  small  as  10 1  H/  (131. 
For  cold  atoms  the  sensitivity  increases  linearly  with  the 
increased  atomic  lifetime  in  the  interaction  zone. 

(3)  Inertial  effects  can  significantly  modify  the  phase  of 
the  GSE.  We  calculate  the  effect  of  a  gravitational  ac¬ 
celeration  along  the  k  direction  on  the  amplitude  of  the 
GSE.  taking  into  account  the  temporal  phase:.  d|.  d;,  dj. 
and  04  of  the  excitation  pulses.  We  find  that  the  phase  of 
the  echo  signal  d,  can  be  expressed  as 

d<  —  d4“  ±  [(d: —  di )  ~  i  kgT$z)  (9) 

where  g  is  the  acceleration  of  gravity  and  +  and  —  refer 


(8) 

I 

to  the  P  and  P'  terms,  respectively.  We  have  ignored 
terms  of  order  kgT}\  and  kgT ij  Du  The  phase  0S  may  be 
measured  by  a  homodyne  technique  using  a  reference 
laser  as  both  the  source  of  the  excitation  pulses  and  as  a 
local  oscillator  for  the  echo  phase  detection.  It  is  clear 
from  Eq.  (9)  that  the  reference  laser  need  be  phase 
coherent  only  over  a  time  of  order  D|.  and  not  over  the 
duration  of  the  complete  echo  sequence  T 42.  Since,  ac¬ 
cording  to  Eq.  (9),  T 21  may  be  made  arbitrarily  short 
without  loss  of  sensitivity  to  g,  a  highly  stabilized  laser  is 
not  required  for  the  measurement.  By  comparison,  the 
gravitational  phase  shift  for  the  SE  is  0,—  d4=±l(o: 
—  di )  —  4g/'42Di].  smaller  than  that  of  the  GSE  by  a 
factor  of  Di/D2~  I / y D 42-  To  achieve  high  precision 
with  the  SE,  a  long-lived  excited  stale  is  required  as  well 
as  a  very  stable  laser  source. 

As  an  example,  we  calculate  the  sensitivity  to  gravita¬ 
tional  acceleration  of  the  GSE  applied  to  a  gas  of  atoms  at 
room  temperature.  Assuming  an  uncertainty  in  the  phase 
measurement  of  <5d,.  the  corresponding  relative  uncertain¬ 
ty  in  a  measurement  of  g  is  Sg/g  =2SoJkgTi2-  For 
k  —  I05  cm-1  and  D:~I0-4  s  (limited  by  the  atoms' 
lifetime  in  the  interaction  volume),  one  finds  Sg/g  —  2<5o,. 
Thus,  for  a  phase  sensitivity  of  i0-'.  one  could  measure 
the  gravitational  force  on  a  gas  of  room-temperature 
atoms  to  a  precision  of  about  0.  lc".  For  a  gas  of  laser 
cooled  atoms,  one  should  be  able  to  achieve  *he  same  sen¬ 
sitivity  5g  as  that  reported  in  Ref.  [7], 
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